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Abstract 

The individual-based models constitute a set of widely implemented tools to analyze the incidence 
of individuals heterogeneities in the spread of an infectious disease. In this work we focus our 
attention on human contacts heterogeneities through two of the main individual-based models: 
mobile agents and complex networks models. We introduce a novel mobile agents model in which 
individuals make displacements with sizes according to a truncated power-law distribution based on 
empirical evidence about human mobility. Besides, we present a procedure to obtain an equivalent 
weighted contact network from the previous mobile agents model, where the weights of the links 
are interpreted as contact probabilities. From the topological analysis of the equivalent contact 
networks we show that small world characteristics are related with truncated power-law distribution 
for agent displacements. Finally, we show the equivalence between both approaches through some 
numerical experiments for the spread of an infectious disease. 
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I. INTRODUCTION 



The spread of a disease in human populations had initially been studied by means of 
deterministic and stochastic compartmental models in which individuals are classified ac- 
cording to their status as susceptible (S), infective (I) or refractory (R) The original 
approach implicitly introduces the assumption of totally mixed population, where all the 
agents can interact each other, that allows a simple mathematical formulation {2, 3|. Sub- 
sequently, the need to incorporate spatial and population heterogeneities has led to some 



reformulations 



The path recently followed to describe the social interactions that lead to infectious 
contacts is the adoption of complex networks as a substrate on which an infectious disease 



spreads . The introduction of complex networks as elements of description has mainly 
provided the incorporation and study of contacts heterogeneities. Thus, in this context, 
each agent interacts with a finite number of neighbors, which in turn varies from agent to 
agent. Many of the efforts have been concentrated on the study of the impact of specific 
network topological features in the spread of an infectious disease jol-lSl. In particular, the 
incidence of degree distribution has been extensively analyzed with respect to occurrence 
probability and total size of an outbreak, especially in some networks with power law degree 
distribution for which has been shown the absence of epidemic threshold for susceptible- 
infective-susceptible dynamic (SIS) 14j-|l7l|. As another example, in 18| the authors have 
found a transition from endemic to epidemic oscillatory behavior for SIRS model in networks 
with Small World property. 

Recent years have yielded notable advances in the study of human displacements driven 
by new possibilities offered by mobile technology systems. In a recent paper [l9[ the authors 
analyze the pattern of movement of a sample of cell phone users finding a kind of truncated 
Levy flight trajectories with a high degree of regularity, contradicting earlier models that 
described the human displacements from an uniform random walk. Althou gh Levy flight 



trajectories have been previously suggested to describe the human mobility [20|, this does 



not necessarily imply that each individual must follow a pure Levy flight displacement 
distribution, as pointed in [19|. It should be noted that Levy flight displacements have been 
extensively observed in nature and, for example, are believed to constitute the optimal food 
search strategy for some animals [2l|, 



221. 
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These findings concerning human mobility allow radical improvements on another widely 
extended modeling tool like mobile agents (MA) models. Its aim is, essentially, seeking to 
replicate human movement and, for example, its impact on the spread of an infectious disease 
transmitted by closeness or direct contact. In [SsJ the authors perform realistic simulations 
of infectious disease transmission in the city of Portland, Oregon, USA. Another MA model 
las been proposed to reproduce the topological features of some networks of acquaintances 



24|. 

On the other hand, while the scientific community has begun to shed light on the impact 
of different topological features of networks in the spread of human infectious diseases, little 
is known about the origins of such features in social contact networks. 

In this paper we consider the spread of a close-contact infectious disease with SIR dynamic 
through a population of mobile agents. With this aim we propose a set of MA models 
considering different movement patterns for the intervening agents. In particular we consider 
two alternatives: one where each agent performs an uniform random walk, and the other 
one that includes recent findings about human mobility, with displacements sizes that follow 
a truncated power law distribution as discussed in [l9[. On the other hand, we propose an 
alternative and, as shown later, equivalent complex network-based approach. For this, we 
show a method to build an equivalent weighted contact network where the weight of each 
link indicates its existence probability. 

After introducing both approaches we focus on drawing a correspondence between the 
displacement pattern of agents and topological characteristics of the associated contacts 
networks. From this analysis we find that the Small World characteristic of the equivalent 
contact networks is related to the agents displacements whose sizes follow a truncated power 
law distribution. Finally, we analyze the impact of periodicity in agent displacements in 
order to reproduce recurrent behavior {25}]. 

The rest of this paper is organized as follows. In Section [III we present the details of our 
MA model as well as the criterions of contact and infection transmission between agents. In 
Section Hm we introduce the procedure to build an equivalent weighted contact network and 
propose alternative quantities to describe its topological characteristics. Section HVl presents 
some numerical examples that show the equivalence between MA models and its associated 
contact networks. Also we analyze the topological characteristics of the equivalent contact 
networks and discuss the link between usual topological features of some social networks 
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and agents mobility pattern. Finally, Sec. |V] summarizes our conclusions and suggestions 
for future works. 



II. MOBILE AGENTS MODEL 

Here we propose a stochastic MA model for the transmission of an infectious disease 
applied on a population of mobile agents, disposed on a square shaped cell of linear size 
L. Initially, agents are ordered on a regular square lattice with distance a = L/(\/N + 1) 
between nearest neighbors. The position for each agent in this initial configuration will be 
named "home" from now. 

In particular, we consider two alternative patterns of motion for individual agents: 1) 
displacements follow truncated power law distribution or 2) uniform random displacements. 
Regarding scenario 1), recent studies from the trajectories of a significant number of cell 



phone users have shown that human mobility has a high degree of regularity [25|. It was 
noted that the distribution of displacements of the analyzed sample can be well approximated 
by a truncated power law [3]: 

P(Ar) oc (Ar + Aro)"^ exp(-Ar/K) (1) 

where /3 = 1.75 ± 0.15, Aro is a characteristic displacement length and k a cutoff value. 

size of the displacements of each 
individual to a time t, namely the radius of gyration r^, also follows a truncated power law: 

P{rg) (X {vg + rj)-^-- expi-rj^r) (2) 



with /3r = 1.65 ±0.15 [l9j. Under these assumptions, they conclude that all agents are char- 
acterized by a single distribution of step size independent of their r^. Then, the distribution 
of displacements for an agent Prg{Ar) can be represented as Prg(Ar) = r~°'F{Ar /vg), where 
F is a scaling function satisfying F{x) ~ for x < 1 and a rapid decreasing for x ^ 1. 
In this scenario, we assign to each agent i = 1, ...,N a radius of gyration r* (t) from Eq. 
with K, = L = 1 and = a. Then, the agents perform displacements with randomly 
chosen direction whose sizes are distributed according to a power law so that its radius of 
gyration r* (t) corresponds to the values previously assigned. 
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Our second alternative for agent displacements is based on the usual simplified repre- 
sentation of individuals as random walkers. As in the previous case each agent chooses a 
random direction, but now the displacement sizes are uniformly distributed within a given 
interval. We show examples from both scenarios in Section IIVI 

A. Displacement probability 

As is well known, some people make more movements per day than other: it is highly 
probable that a postman makes more displacement a day than an office worker. In order to 
incorporate these heterogeneities in our model we have assigned each agent a displacement 
probability per step through a predetermined probability distribution. We focused on two 
usual distributions: i) Poissonian and ii) Power-Law type distribution Pi{x) ~ x''' , where 
X represents the number of movements per day made by agent i. For example, if Pj = 1, 
the agent i will make a displacement for every step of the temporal evolution. As regards 
the time evolution, we consider a discrete evolution taking 1 hour as our time unit, in other 
words, every step of our algorithm represents an evolution over time of 1 hour. Then, the 
maximum number of movements per day for a given agent will be x = 24. Finally, at the 
end of each day the agents return home with constant return probability Pret = 0.9, in 
order to consider some trips more prolonged in time too. 

In order to keep constant the number of agents in the cell we establish periodic boundary 
conditions for the movements but not for interactions, i.e, we do not consider periodic 
boundary conditions to calculate distance between agents. This means that when an agent 
leaves the cell through one of its sides, he will be re-entering by the opposite one, however 
two agents located in opposite sides do not interact. 

B. Periodic movement pattern 

As human beings we have acquired habits which lead us to reiterate our movements 
with some frequency, for example: some of us take the kids to school in the morning, 
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then go to the office, less frequently go to gym and finally return home. It is possible 
that sometimes we perform specific non-routine activities but essentially our trajectories 
are predictable [25]. Therefore, in order to reproduce this behavior, our model incorporates 
an independent periodic motion pattern for each agent. This means that each agent i = 
1, ...,N periodically repeats a given sequence of qi positions r'j^, rg, r^. consistent with 
previously established scenarios. This is achieved by assigning to each agent i = 1,...,N a 
period tj. In our particular case, the set {ti}i=i,...,Ar is uniformly distributed in [l,Tp], where 
Tp is the maximum allowed period. 

C. Infectious disease transmission 

Now that we have detailed the features of our mobile agents model, we will focus on 
the transmission process of an infectious disease. We study the spread of a general and 
schematic disease corresponding to Susceptible-Infectious-Recovered (SIR) dynamic. The 
transmission occurs by close proximity through direct or indirect contact from infected to 
susceptible agents, as for infiuenza, measles or chickenpox, although we will not focus on 
the specifics features of any given disease. Here we want to present a general framework 
to eventually incorporate the specificities of a given infectious disease transmitted by close 
contact. 

At this stage it is necessary to establish a proximity criterion between agents that leads 
to the transmission of an infectious disease. In this model we consider that physical contact 
is not the only way to produce contagion, but just enough that a susceptible agent is near an 
infected one, sharing a given physical space. This assumption is based on the transmission 
processes of some infectious diseases, like infiuenza, which can be transmitted by means of 
Fliigge's droplets emitted by a sneeze of an infected individual. 

We start the process with a single infected individual located at the center of the square 
cell. After updating the position of all agents to perform one evolution step, a contact 
radius is established around its positions. Thus we say that every pair of agents at a 
distance less than Vc are in contact or linked. Then, the probability that a susceptible agent 
i becomes infected during an evolution step is given hj Pi{zi) = 1 — {1 — PcY\ where Zi is 
the number of infected individuals found at distance less than from agent i and pc is the 
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infection probability per contact. Each infected agent remains in that state for Tj steps of 
time evolution, after which it becomes recovered (R) and then immune to disease. It should 
be noted that the updating process is performed synchronously. 



III. EQUIVALENT COMPLEX CONTACT NETWORKS 

So far we have described a MA model where the disease transmission is associated with 
the movement pattern followed by agents. This may help us to understand how spatial 
displacements of people affects the spreading of an infectious disease. 

Another alternative path followed in the literature to incorporate spatial heterogeneities 
in epidemiological modeling is to introduce a complex network as a substrate on which 
a disease spreads. Thus, the network nodes represent individuals, while the edges their 
contacts or relationships. There is abundant literature about infectious disease spreading on 



contact networks with diverse topological features [9l-ll7l]. However, there are few cases in 
which it has been shown the correspondence between the studied system and the proposed 
contacts network. 

In this section we will show how an equivalent weighted and undirected contact network 
can be drawn from the MA model described in previous section. 



A. Network construction. Time window 



As stated in our proximity criterion in Section Hi C\ we say that two agents are in contact 
when they are at distance less than re- In this case, those nodes, say i and j, will be joined 
by a link lij whose weight Wij ^ 0. The set of nodes representing agents, together with their 
links are what we call equivalent contact network. For the purpose of its construction we 
need to register contacts between agents during a time interval or time window Ty/. If all 
agents have a recurrent behavior or, still better, a periodic motion pattern, a given agent i 
repeats its movements with period tj, visiting the same sites over and over again. 

In this case Tw must be long enough to allow all agents to complete their movement 
patterns, i.e., 

Tw > max (tj) (3) 

i=l,...,A'' 

However, in Section |IV] we will show through some examples that the presumption of recur- 
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rent behavior in agents movements is not a necessary condition for the construction of an 
equivalent contact network. It should be noted that for non-recurrent displacement patterns 
it becomes necessary to consider the infectious period Tj as an additional constraint on Tw 
value, resulting Tw > Ti. 



B. Probabilistic interpretation of weights 

Following our consideration about network construction, it is clear that some contacts 
will be more frequent than others, so that the resultant contact network will be weighted. 
We know the identities of all contacts for every agent i over Tw and also the num- 
ber of times they meet each other. Then, link weight Wij between nodes i and j can 
be given by the number of contacts per unit time between these nodes during Tw, i.e., 
Wij = (number of contacts between i and j during Tw)/Tw- (There exist other factors that 
contribute to the heterogeneity of contacts among agents, like the distinction by age group, 
race or previous history of disease, neglected in this study with respect to contacts fre- 
quency.) However, here we give a probabilistic interpretation of weights. If we define nij as 
the total number of contacts between nodes i and j during Tw, then, the elements Wij of 
weight matrix W are given by: 

During each time evolution step nodes establish contact only with a subset of their 
neighbors, chosen with probabilities given in Eq. HJ Thus the actual pattern of contacts 
makes stochastic changes over time in order to reproduce the behavior of MA model. As 
a result of the described approach we obtain a matrix W that contains the fundamental 
features of the equivalent contact network, such as its undirected character evidenced in Eq. 
il 



C. Infection process in the equivalent contact network 

The infection spreading process in terms of the equivalent contact network is as follows: 

1. We start the spreading process with one infected agent placed at the center of the cell, 
as we do for MA model. Each agent i = I, ...,N has an assigned counter tc{i) that 
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defines its state as usual for SIR infectious diseases: 

tc{i) = agent z is S 
1 < < Ti agent i is I (5) 
tc{i) > Ti agent i is R 

2. Each susceptible node makes contact witli a given set from its neigliborhood according 
to contact probabilities established by matrix W. 

3. Now we apply the same infection criteria as for MA model: Suppose that a given node 
i makes contact with n* nodes from its neighborhood, Zi of which are infected. Then, 
i will become infected with probability Pi{zi) = 1 — (1 — Pc)^"- 

4. We perform the contagion criteria synchronously over all S agents. 

5. The counters of all infected agents are actualized: tc{i) = tc{i) + 1 if agent i is in 
infected state (I). 

6. The states of the infected agents are actualized so that each infected agent i remains 
in that state if tc{i) < Ti or it makes a transition to recovered state (R) if tc{i) > Ti. 
This completes one evolution step. 

Figure [1] shows an example of an equivalent contact network obtained for N = 10^ mobile 
agents following a truncated power-law displacements distribution over a time window of 
Tw = 7 days. 

Now we present a reformulation of two classical topological quantities, i.e., average path 
length {£) and clustering coefficient (C), in order to characterize this weighted contact net- 
works. 

D. Average path length and clustering coefficient in weighted contact network 

One of the quantities that characterize the topology of complex networks is the average 

Bith length. It can be defined as the average minimum distance between any pair of nodes 
, Q] and, from now, we call it the topological average path length {£). When we refer 
to distance in this case we mean the minimum number of links that must be traversed to 
communicate a given pair of nodes. There exist several algorithms in the literature to find the 
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FIG. 1: (Color online) Equivalent contact network obtained from a MA model with N = 10^ agents 
following a truncated power-law displacements distribution. Each node has been drawn with a size 
proportional to its degree and with a color representing its state: susceptible (green), infected (red) 
or refractory (blue). 

shortest-path in non-weighted networks, some examples are Dijkstra's algorithm. Breadth- 
first search and Floyd- Warshall algorithm. In addition, some of them can be generalized for 
weighted networks. 

Generally, for ordinary weighted networks, each link has an associated weight corre- 
sponding, for example, to its ability to carry information or, as in the problem studied here, 
spreading a disease. According to our probabilistic interpretation, each weight indicates the 
probability of link existence. Then, a given path leading from node i to j through the links 

km, lnun2, ^na.ng, ••• , ^n„,j wiU haVC existCUCe probability Pi^j = Win^Wnin2U!n2n3---Wn,^j- 

This is why it is necessary an additional definition of average path length for this probabilis- 
tic interpretation of weights. The alternative definition should include both the distance d^j 
between nodes i and j through path k as the existence probability of that path. Then, we 
call it the probabilistic average path length (ip). 

The probability Pjl^j can be interpreted as a measure of the persistence of the path k 
that connects node i to j, so that the probabilistic shortest path is the one that maximizes 
the probability Pjl^j- In this sense, we say that dij = rf^- is the minimum distance between 
i and j if: 

Pl'-^,= max (6) 

k£{all paths ■' 
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where k runs through all paths connecting nodes i and j and we call L to the probabilistic 
shortest path that connect them. Finally, the probabilistic average path length (ip) will be 
the average of d^j between all possible pair of nodes in the network: 

1 ^ 

Another quantity used to characterize network topology is the average clustering coef- 
ficient C. It measures the likelihood that two neighbors of a node i are associated among 



themselves 



71]. The local version of clustering coefficient c, can be defined for every node 



i as the fraction of connected neighbors of i, so that the average clustering coefficient is 

i 

If we call n{i) the set of ki neighbors of i and we define Wi = Wij, then the local 

clustering coefficient can be written as: 



The definition of Eq. M applies to both weighted and non- weighted networks since when 
Wij = 1 Wi,j = 1, usual definition for non-weighted networks is recovered. 



IV. NUMERICAL SIMULATIONS 

We consider a mobile agent model with A^ = 10"^ agents within a square cell of side L = 1 
(in adimensional units). We fix the contact radius Vc equal to half the distance between 
nearest neighbors in the initial square lattice configuration, i.e. rc = a/2 = L/ {2{\/N + 1)). 
Actually, the value of Tc will depend on the disease characteristics but, as already mentioned, 
we are not concerned here with any particular disease. 

The displacement probability for each agent was assigned from a Poisson distribution 
P{x) = ^e-^ with A = 2.1. 

Each agent performs a periodic motion with step size according to the truncated power 
law distribution described in Sec. Ull We will refer to this model as periodic truncated power 
law displacements (PTPL). The periods of movement {ti}i=i,...,7v are chosen with uniform 
probability from [1, 7], i.e. from one day to a week. 
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FIG. 2: (Color online) Contact network topology characterization for NPUR (red squares), PTPL 
(blue circles) and NPTPL (green triangles) variants of our MA model, a) Average probabilistic 
path length versus time window size T^. b) Clustering coefficient for different Tw. Each point 
correspond to average over 1000 contact networks derived from independent realizations of each 
variant of MA model. 

A. Time window election and topological properties of contact networks 

Now we will construct an equivalent contact network from PTPL model. For this we 
need to determine an optimal time window size T^- Common sense dictates that the larger 
T\Y the better will be the description given by the equivalent contact network. Clearly there 
is a practical compromise, so it is necessary to establish a lower bound for Tw so that the 
topological characteristics of the equivalent contacts network are not substantially modified. 

As can be seen in Fig. [2] for PTPL, both ip and C approach their asymptotic values 
for Tw > 7, consistent with the maximum period of our mobile agents. It should be noted 
that £ and £p are approximately equal for the studied cases and therefore we do not show 
both results. The equivalent contact network obtained for Tw = 7 has average degree 
(k) = 8.3 ± 0.6 and Small World topological characteristics with i = ip = 7.12 ± 0.38 and 
C = 0.374 ±0.004. 

It is important to note here that the Small- World features arise naturally from the impo- 
sition of a movement pattern for agents, with displacements sizes according to a truncated 
power law distribution. 

On the other side, the incorporation of periodicity in agent movements limits the spatial 
locations available to each agent, reducing the time window necessary to obtain a correct 
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FIG. 3: (Color online) Average degree {k) in terms of time window size T^y for equivalent contact 
networks obtained from NPUR (red squares), PTPL (blue circles) and NPTPL (green triangles) 
mobile agents displacement models. Again, each point in the graphic corresponds to average over 
1000 independent realizations of each variant of MA model. 

sampling of possible contacts between them. To verify this last assertion we repeat the 
analysis above but now for two alternative scenarios: i) agents following non-periodic trun- 
cated power law (NPTPL) displacements, and ii) agents carrying out non-periodic uniformly 
distributed random (NPUR) displacements. 

As shown in Fig. [21 NPUR model has slower convergence to the asymptotic regime com- 
pared with PTPL model, in particular decaying as a power law. This behavior is consistent 
with the linear dependence between average degree and Tw observed for NPUR model in 
Fig. [3l In contrast, when agents follow a periodic movement pattern as in PTPL, the curve 
representing the evolution of {k){Ty/) also tends to an asymptotic regime. The NPTPL 
model shows an intermediate behavior between the previous ones. Clearly, when Ty/ = 1 is 
difficult to distinguish between a model with agents making uniform random displacements 
(NPUR) and another whose agents follow a truncated power-law step size distribution with- 
out any periodicity (NPTPL). By increasing Ty/ the contact network associated to NPTPL 
model begins to show the Small World characteristics also observed for PTPL displacements 
model (see Fig. [2]). 

On the other hand, setting Ty/ = 15 for NPUR model we obtain an equivalent contact 
network with i = ip = 3.1±0.1 and C = 0.014 ±0.001, characteristics consistent with 
Random Networks. 
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FIG. 4: (Color online) Degree distribution for contact networks derived from NPUR (red squares), 
PTPL (blue circles) and NPTPL (green triangles) mobile agents models. In all cases we consider 
T\Y = 7 and we average over 1000 independent realizations of each variant of our MA model. 

Despite the foregoing, in Section |V] we show that the disease spreading dynamic for the 
NPUR model can also be reproduced with very good matching through its equivalent contact 
network. 

Figure H] shows the degree distributions for PTPL, NPTPL and NPUR models, all for 
Tw = 7 days. For PTPL the degree distribution initially behaves like a Poisson distribution 
but then decays more slowly, following a power-law with exponent r = —3.2 for 10 < A; < 35. 
This heavy-tailed degree distribution has been previously reported for urban social contact 



networks 



In contrast, degree distribution for NPUR shows an irregular behavior resulting in an 
increased frequency of high-degree agents. NPTPL resembles the behavior of NPUR for 
1 < < 5 but for k > 5 it behaves like PTPL. 

B. Infectious disease spreading. Comparison between Mobile-agents model and 
its equivalent contact- network. 



We will now analyze the spread of an infectious disease with SIR dynamic in a population 
of PTPL mobile agents, in order to compare its results with those obtained from the equiv- 
alent contact network approach. We select a fixed infectious period Tj = 4 for all agents and 
a time window T^r = 7, that satisfies T\y > maxj=i Ar(tj) and T]v > ti. Figure E] shows 
the time evolution of total infected agents for PTPL mobile agents model and its equivalent 
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FIG. 5: Dynamical comparison between the MA model and its equivalent contact network approach 
for the spread of an infectious disease with SIR dynamic. The time evolution of total infected 
individuals for PTPL mobile agents model is represented by the solid line in the figure, while the 
squares correspond to the results for its equivalent contact network. All the results were obtained 
from the average over 10000 independent realizations of the disease spreading process on a single 
spatial realization of the PTPL mobile agents model and its equivalent contact network. The inset 
shows the degree distribution for the contact network with average degree {k) = 8.905, obtained 
from the imposition of a poissonian displacement probability for the agents in MA model. 

contact network. It can be seen that both approaches show a very good agreement. 

Equivalently, in Figure |6] we compare the total number of infected agents for NPUR 
mobile agents model with results for equivalent contact networks corresponding to different 
time window sizes. The good agreement observed for T^y > 6 in Fig. [H] can be understood 
remembering the constraint Tw > Tj and that, in the context of NPUR model, each mobile 
agent makes an uniform random walk, thereby increasing its number of contacts with time, 
thus improving the sampling performance. 

So far we have arbitrarily assumed that the displacement probability, i.e. the probability 
that an agent makes a displacement, follows a Poisson distribution. Now we show that the 
observed agreement between both approaches is not constrained to this condition by impos- 
ing a power-law displacement distribution P{x) ~ (x — xo)~'^ ex.p{—x/ko) with exponential 
cutoff. All parameters were selected so that the average degree is close to that obtained 
for the example shown in Fig. [51 The equivalent contact network also shows Small- World 
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FIG. 6: Comparison of the total infected agents evolution for a set of equivalent contacts networks. 
The sucession of equivalent contact networks has been obtained by increasing the sampling time 
Tw from 1 to 7 days over a single spatial realization of the NPUR mobile agents model. The solid 
line shows the time evolution of infected individuals for NPUR mobile agents model. Again, all 
the results correspond to averages over 10000 independent realizations of the disease spreading 
process. 

topological characteristics with £ ~ £p = 6.977, C = 0.383 and average degree {k) = 8.324. 
Again, we show in Fig. [7] the good agreement between PTPL mobile agents model, in 
this case with power-law displacement probability, and its equivalent contact network, for 
temporal evolution of infected agents. 



V. SUMMARY AND CONCLUSIONS 



In this paper we introduce a novel MA model in which agents make displacements accord- 
ing to a particular truncated power-law distribution consistent with recent findings related 



to human mobility |l9|, |20| . We have developed a procedure to build an equivalent weighted 



contact network by recording contacts through a proximity criteria between mobile agents 
during a given time window. The weights are interpreted here as contact probabilities and 
this fact provides stochasticity to the contact network. Given the probabilistic nature of the 
weights we have introduced alternative definitions for the topological network descriptors. 
In particular we have defined a probabilistic average path length which results equivalent 
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FIG. 7: Time evolution of the infected individuals for PTPL mobile agents model with power law 
displacement probability (solid line) and its equivalent contact network (squares). The inset in 
the figure shows the degree distribution of the equivalent contact network, which now resembles a 
power-law behavior. 

to the usual topological average path length in the studied cases. We want to note that the 
analyzed contact networks are not necessarily friendship networks or coworkers networks but 
just are networks of people who have shared the same spatial region for a period of time. 
However, contact networks can also provide information about friend or work relationships 
among individuals through the weights of the links that bind them. 

We have analyzed three different agents movement patterns for which we have studied 
the topological features of its equivalent contact networks. Our results provide evidence that 
relates the truncated power-law displacement pattern of agents with Small World character- 
istic of its equivalent social contact networks. The former correspondence can be explained 
by considering that this kind of movement pattern increases the encounters probability 
among agents whose homes are close with respect to those who are far away. In this way, 
there will be tightly interconnected sets of agents, contributing to increase the average clus- 
tering coefficient, which are linked to others by sporadical long displacements that act as 
shortcuts. We are currently addressing the study of community structure formation from 
different displacement patterns and its impact on spreading dynamic in complex networks. 
Another important factor to consider, related to the previous one, is the distribution of 
homes in the cell. In this paper we have considered a regular and uniform distribution of 
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homes in order to reduce the number of parameters involved in the problem. However, this 
is an appropriate assumption only for some urban areas, so it is necessary to extend our MA 
models to study the impact of demographic heterogeneities. 

The contacts networks obtained from PTPL model with regular reincidence in agent 
positions reached its asymptotic regime, with respect to topological variables, faster than 
the other two displacement models without reincidence. This assumption of regularity in 



human mobility has recently received empirical support [25| , allowing us to optimize the time 
of sampling that leads to equivalent contact networks. However, we have shown through an 
example for the spread of a SIR infectious disease with finite infectious period, that a MA 
model with non-recurrent movement pattern as NPUR model shows good agreement with 
the results of its equivalent contact network. 

Finally, we have shown the equivalence between mobile agents and contact networks 
models for the spreading of a particular kind of infectious diseases with SIR dynamic and 
transmitted by proximity. 
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